10/523124 



Attorney Docket No. 2003UR01 1 



m 



BTOI 




d PCT/PTO 2 6 JAM 2M5i 



METHOD TO CONVERT SEISMIC TRACES INTO PETROPHYSICAL 

PROPERTY LOGS 

[0001] This application is the National Stage of International Application No. 
PCT/US04/04313, filed February 13, 2004, which claims the benefit of U.S. 
Provisional Patent Application No. 60/458,093, filed March 27, 2003. 



[0002] This invention relates generally to the field of reservoir characterization 
using petrophysical properties. More specifically, this invention relates to deriving 
petrophysical property data with seismic data as constraints. 



[0003] Reservoir characterization is a process of using well log and seismic data 
to map reservoir geometry, net-to-gross ratio, porosity, water saturation, permeability, 
and other petrophysical properties. Each seismic response trace is the result of a 
seismic source wavelet convoluted with a spatial earth reflectivity series, if linear, and 
is not a direct indicator of rock and fluid properties. Well logs with much higher 
vertical resolution than seismic traces, such as Gamma Ray, Density, Porosity, Sonic 
velocity. Resistivity etc. measured from borehole can be direct indicators of rock and 
fluid properties after some borehole environmental corrections. 

[0004] The number of wells drilled in a reservoir is always limited because of well 
cost and the thickness of many reservoir formations is often below seismic resolution. 
If wells were drilled in each location of the seismic trace (or each location of 
geophone), then logs would be obtained from these wells and the reservoir 
characterization would utilize the well logs to map reservoir geometry, net-to-gross 
ratio, porosity, water saturation, permeability, and other petrophysical properties. For 
economical reasons (or cost-effectiveness), there are rarely enough wells to properly 
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map the details of reservoirs. Therefore, Petrophysical properties are inferred form 
the seismic data. 

[0005] Until recently, seismic attributes were extracted and calibrated with 
available well control data, and petrophysical properties were 
interpolated/extrapolated between and beyond sparse well control, A set of seismic 
attributes related to time, amplitude and shape, including trace-to-trace time tracking 
cross correlation coefficient, integrated absolute amplitude, and area ratio of fatness, 
were calculated and extracted from reflection seismic trace(s). Literval attributes such 
as, the average of all positive peak amplitudes and integrated absolute amplitude are 
calibrated against petrophysical properties from well logs. Rock property maps can be 
generated utilizing these attributes at an interval average. Though many attributes can 
be generated from a seismic traces, how to classify and calibrate attributes with well 
data for best quantitative reservoir characterization and uncertainty analysis are 
continual subjects of research. 

[0006] U.S. Patent No. 5,444,619 (Hoskins et al.) uses a "Artificial Neural 
Networks (ANN's)" method to estimate the relationship of reservoir properties and 
seismic data.. However, the confidence level of such prediction using statistical 
methods are often questioned and challenged due to lack of basic physical 
relationships between a seismic trace and petrophysical properties indicated by a set of 
logs. Another method for determining reservoir properties from seismic data is 
disclosed in U.S. Patent No. 5,835,883 (Neff et al.) This method requires constructing 
pseudo well logs and synthetic seismic traces of the seismic and well log data to 
produce an estimate of the petrophysical properties. The model includes a matching 
scheme wherein trend curves of petrophysical properties are selected based on the 
match to the data. This process requires operator selection and selecting the wrong 
trend curve can lead to errors. 

[00071 Accordingly, there is a need for an innovative method to directly convert 
seismic data to petrophysical property data. The methods preferably includes 
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converting three-dimensional seismic data to petrophysical property cube(s) in either 
time or depth domain with enhanced vertical resolution. Embodiments of this 
invention satisfy these needs. 

SUMMARY 

[0008] A method to convert seismic traces to petrophysical properties is disclosed. 
The method comprises (a) deriving a Combined Log-Seismic Response (CLSR) filter 
from at least one petrophysical log from a well and at least one seismic trace near the 
well and (b) convolving the CLSR filter to convert each seismic trace to a "pseudo" 
petrophysical log. 

[0009] A second embodiment of the method converts a three-dimensional seismic 
volxmie cube of an area to a three-dimensional petrophysical property volume cube. 
The method comprises (a) providing a three dimensional seismic volume cube from a 
seismic survey of the area, (b) providing at least one petrophysical log for the well in 
the area, (c) deriving at combined log seismic response filter from a log and at least 
one seismic trace near the well, (d) converting seismic wiggles traces to petrophysical 
properties by convolving the CLSR filter, (f) outputting the petrophysical properties as 
a three-dimensional property volume cube. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0010] Fig. 1 is a flow chart illustration of one embodiment of a method to 
convert a seismic wiggles to petrophysical properties; 

[0011] Fig. 2 is a flow chart of one embodiment of a method to convert a three- 
dimensional seismic survey and at least one petrophysical well log in the seismic 
survey to a petrophysical property cube; 
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[0012] Fig. 3 is a graphical flow chart of an embodiment of the invention; 

[0013] Fig. 4(A) is a high-resolution well log response function in time domain; 

[0014] Fig. 4(B) is the frequency spectrum of the well log response function in 
Fig. 4(a); 

[0015] Fig. 4(C) is a low-resolution seismic wavelet; 

[0016] Fig. 4(D) is the frequency spectrum corresponding to the seismic wavelet 
in Fig. 4(C); 

[0017] Fig. 4(E) is a high-resolution petrophysical log from a well; 

[0018] Fig. 4(F) is the frequency spectrum of the high frequency well log from 



[0019] Fig. 4(G) is a low-resolution seismic trace near the borehole from which 
the well log in Fig. 4(E) was run; 

[0020] Fig. 4(F) is the frequency spectrum of the low frequency seismic trace 
corresponding to the seismic trace in Fig. 4(G); 

[0021] Fig. 4(1) is the Combined Log-Seismic Response (CLSR) from the seismic 
trace of Fig. 4(G) and the well log of Fig. 4(E); 

[0022] Fig, 4(J) is the frequency spectrum of the Combined Log-Seismic 
Response (CLSR) of Fig. 4(1); 
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[0023] Fig. 4(K) is the predicted high-resolution petrophysical log obtained by 
convolving the CLSR of Fig. 4(1) with the low-resolution seismic trace of Fig. 4(G); 

[0024] Fig. 4(L) is the frequency spectrum of the predicted high frequency log of 
Fig. 4(K). 

DETAILED DESCRIPTION 

[0025] In the following detailed description, the invention will be described in 
connection with its preferred embodiments. However, to the extent that the following 
description is specific to a particular embodiment or a particular use of the invention, 
this is intended to be illustrative only. Accordingly, the invention is not limited to the 
specific embodiments described below, but rather, the invention includes all 
alternatives, modifications, and equivalents falling within the true scope of the 
appended claims 

[0026] The present invention is an iimovative method to convert a seismic trace to 
a petrophysical property log property and generate petrophysical property data 
including for example, three-dimensional cubes in either time or depth domain with 
enhanced vertical resolution up to log resolution. The method is based on rigorous 
wave physics and petrophysics principles. 

[0027] The high resolution petrophysical logs as traces corresponding to each 
surface seismic trace can be derived by convolution of a Combined Log-Seismic 
Response (CLSR) filter with each low resolution seismic trace. This allows 
generation of two-dimensional slices and three-dimensional volumes of petrophysical 
property data. As shown in figure 1, one embodiment of the invention includes two 
steps. First, a Combined Log-Seismic Response ("CLSR") filter is derived 1, The 
CLSR filter is derived from a petrophysical log and at least one seismic trace near the 
well. In the second step, the CLSR filter from the first step is convolved with each 
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seismic trace to convert each seismic wiggle trace to petrophysical properties 3. 
Converting a cube of three-dimensional seismic volume traces will generate three- 
dimensional petrophysical property volume cubes. 

[0028] In the seismic data, the high-frequency reflectivity information are usually 
filtered out due to the nature of the low-frequency seismic wavelet, and the high- 
frequency information cannot be brought back through conventional seismic 
processing. However, when both a high-resolution log and a low-resolution seismic 
trace are available at the same spatial location, a relationship can be determined 
between them. Then, the relationship, described below, may be used to extrapolate 
(predict) the high-resolution logs at a neighboring location where only the low- 
frequency seismic data are available. The extrapolation requires an assumption that 
the relationship is mathematically stable and does not vary spatially within the survey 
area or reservoir unit. 

[0029] Following one embodiment of the convolution model, in time domain, a 
linear response of a low frequency seismic response lseis(t) and a high frequency well 
log response hiog(t) may be expressed as the convolution of the low frequency seismic 
response filter Wi(t) and high frequency logging tool response filter Wh(t) with the 
petrophysical property series at the same spatial location as expressed in Equation (1) 
and (2). 

[0030] A low frequency seismic response lseis(t) and a high frequency well log 
response hiog(t) may be expressed as the convolution of the low frequency seismic 
response filter wi(t) and high frequency logging tool response filter Wh(t)^ The 
petrophysical property series at the same spatial location may be expressed as 
Equations (1) and (2): 



lseis(t) = r(t) * Wi(t) 



(1) 



hiog(0=x(t) * wUO 



(2) 
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where r(t) is the reflectivity of the formation at or near wellbore and x(t) is a 
petrophysical property log converted in time series for the well. Examples of 
petrophysical property logs include but are not limited to gamma ray \02rGR(t), 
resistivity log-iffO, porosity log-PHIft), density log-RHOB(t), or sonic log-DT(t) etc. 
from the well. Furthermore, "*'* stands for convolution, Wh(t) and are high- 
resolution log response function and low-resolution seismic wavelets, respectively. 

[0031] hiog(t), lseis(t\ ^h(Of and wi(t) are respectively waveforms of logs and 
seismic traces and can be expressed as "relative" amplitude as a function of time 
(typically in miUiseconds), and reflectivity r(t) is imitless. Typically the logs have the 
following units, GR is in API, Rt is in OHMM, PHI is in P.U., RHOB is in g/cm^ DT 
is in microsecond/ft. 

[0032] In taking the ratio of Equations (1) and (2), the relationship between the 
log and seismic responses maybe expressed as Equation (3): 

hiog(t)^lseis(t) *c(t) (3) 

where c(t) is the Combined Log-seismic Response (CLSR) filter, and may be 
expressed as c{cj) = Wf^{6})/Wi{o}) in frequency domain. piog(t) is the time domain 
response of a petrophysical log x(t) corresponding to the reflectivity series r(t) that is 
derived from density log-ROHB(t) and sonic log-DT(t), which can be expressed in 
frequency domain asP^^^{co)= X{co)/R{co). When solving for the high resolution 

acoustic reflectivity log, x(t) equals r(t), and then, piogft) or goes to one. In the above, 
c(t) is represents a "seismic-petrophysics covariance" filter which combines the low- 
resolution seismic response with the high-resolution log response, namely, the 
Combined Log-Seismic Response (CLSR) filter. 

[0033] In theory, c(t) can be solved by deterministic deconvolution (for example, 
c(t) = h^^ (t)* (r)"' . A stable c(t) from real data with a finite length may only be 

inverted if seismic response filter wift) has broad bandwidth and is inclusive of the 
logging tool response filter Wh(t). There are methods that may be employed to 
broaden the wi(t) bandwidth and provide a stable c(t). For example, a pre- whitening 
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technique includes superimposing an impulse response curve to a filter design but not 
to final data. (Sheriff and Geldart, 1995, Exploration Seismology, p. 295-299). 

[0034] In order to construct a stable c(t)y its denominator term wi(t) should be of 
broad bandwidth inclusive of the nimierator term Wh(t). However, wi(t) usually has 
little or no energy at high fi*equencies. Consequently, it would be unstable to 
determine c(t) in fi-equency domain by spectral division. Though pre-whitening has 
been effective in treating such deconvolution (Yilmaz, 1987, p. 103), it will introduce 
high-fi:equency noise harmfiil to the extrapolation and images, A practical method is 
to solve the deconvolution as a constrained inversion in time or depth domain where 
amplitude and phase information are treated simultaneously. Then the instability 
problem in deconvolution will become a non-uniqueness problem in the inversion. 

[0035] The process in time domain can be carried out in two steps. In the first 
step, assume that over an area, the low-resolution seismic data were acquired at many 
locations: { IseisfOp 7="^ ...k...}, and at the kth location, a well is drilled and high- 
resolution log response, say lseb(t)k was taken near the kth seismic trace location. In 
the first step, the CLSR filter c(t) can be determined firom: 

C(t)^ hiog(t)k * lseis(t)k' (4) 

since both IseisftJk and hiog(t)k are known, the above is a deterministic process. 

[0036] In the second step, assuming that c(t) is invariant over the seismic survey 
area or a sub-area, the high-resolution petrophysical property log responses h 'log (t)j 
firom seismic trace j for all seismic traces f lseis(Op J^^^f ...k...} within the survey 
area, can be extrapolated (or predicted) using Equation 5. 

h %g(t)j = IseUOj * C(t) (5) 

[0037] This is again a deterministic process because both terms in the right hand 
side of (5) are known. When the logging tool responses Amotion Wh(t) run in the wells 
in a seismic survey area, and the seismic wavelet wi(t) for the survey (such as an 
airgun for marine seismic survey) are invariant throughout the survey, the CLSR filter 
c(t) may also be invariant to satisfy the assumption. Otherwise, hiogft) the high- 
resolution log and lseis(0 the low-resolution seismic trace will not be self consistent. 
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Self consistent is when both the logging tool responses function Wh(t) run in the wells 
in a seismic survey area and the seismic wavelet wi(t) for the survey satisfy the 
invariant assxmiption. If the c(t) also satisfies the invariant assumption over the 
seismic survey area, then the high fi-equency well log response hiogft) and the low 
frequency seismic response Iseisft) may be.called self (or intemally) consistent. 

[0038] As shown in Figure 1, the second step requires convolution of the CLSR 
filter 3. There are three factors in this convolution. The three factors are seismic trace 
IseisfOj reflectivity rft), and wavelet, wi(t) as shown in convolution Equations (1) and 
(2). 

[0039] The CLSR filter c(t) may be solved by deconvolution as shown by 
Equation (4), and the conventional deconvolution has been focused on statistical 
estimation of two out of the three factors (or terms) involved in a convolution. (See, 
for example, Robinson and Treitel, 1980; Ziolkowski, A., 1984; Yibnaz, 1987;. 
Similar to Extrapolation by Deterministic Deconvolution (EDD) (Zhou et al., 2001), 
the approach here can yield more robust results. The convolution process is solved as 
a constrained inversion and is described below. 

[0040] The convolution Equation (3) can be represented in matrix- vector forms: 

Gl = Lc ^ h, (6) 

where I = (l(h In h, ..•>)^, c = (co, cj, C2, ...fy h = (ho, hu ...f are coefficient 
vectors of /(t), c(t), and A(t), respectively; the two matrices G and L consist of 
successions of down-shifled columns of c and /, respectively. C is solved in Equation 
7, A is a high-resolution log that was measured in the wellbore, and / is a low- 
resolution seismic trace from the seismic data. 

[0041] Applying Equation (4), the first step is to derive the CLSR filter, that is to 
invert c from the kth location with well logs with Equation (7): 

LkC = hk. (7) 

The singular value decomposition of Lk may be expressed as: 
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Lk = UAV". (8) 

where U and Fare unitary matrices, and A is the diagonal matrix containing singular 
values {Si}, A unitary matrix has its Hermitian conjugate equal to its inverse, and a 
diagonal matrix has the form of a,y=i for all i=j, and Uij^O for i^. Then the Combined 
Log-Seismic Response (CLSR) filter is 

c = VA' U^hk, (9) 

[0042] In one embodiment, the petrophysical property logs are predicted using 
Equation (5). The prediction of the high-resolution petrophysical logs at location y is 
determined by solving for Equations (10), (1 1) and (12): 

h/^AA^b (10) 

where 

A= Lj V (11) 

ft= U'h, (12) 

[0043] A generalized inversion means eliminating some of the smallest singular 
values that cause instability, though solutions of different cutoff levels of the singular 
values may fit data equally well. (Aki and Richards, 1980). It is computationally 
efficient to construct a matrix D whose element at the m-th row and n-th colimin has 
the form 



dmn= anmbn/Sn (13) 

where anm is an element of A, and is an element of A. A solution of At/' with up to p 
singular values is just a row-summation of the first p columns of the matrix D. 

[0044] The instability problem in the convolution can be treated as the non- 
uniqueness problem in inversion, as manifested by the multiple solutions offered by 
the matrix D. Though the solution is non-unique, it is certainly stable. One method to 
choose a solution is to select the cutoff level of the singular value based on three 
constraints. First, the solution should give the best prediction at the kth location where 
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the high-resolution log is given. Second, a low-pass filtered version of the prediction 
should correlate well with the low-resolution data everywhere. Finally, the ratio 
between magnitudes of the prediction and the low-resolution data may be chosen to be 
close to the ratio at the ktk location with known high-resolution data. 

Example 

[0045] Figure 2 illustrates an embodiment of the method of producing three- 
dimensional petrophysical volume cubes using at least one well log and a three- 
dimensional seismic volume. A three-dimensional seismic survey of the area is 
obtained 20. A well is drilled in the area of the seismic area 21. Alternatively, an 
existing well or combination of wells located in the seismic area may be utilized. At 
least one petrophysical log is obtained from at least one well 22. A CLSR filter is 
derived from a petrophysical log and at least one seismic trace substantially near the 
well 23. The seismic wiggle traces from the seismic survey are then converted to 
petrophysical properties by convolving the CLSR filter 24. The petrophysical 
properties may then be outputted as a three-dimensional property volume cube 25. 

[0046] Figure 3 illustrates one embodiment of the conceptual process of the three- 
dimensional petrophysical volume cubes using well logs and three-dimensional 
seismic volume. Within a three-dimensional seismic survey cube 31, a well 35 is 
drilled and high frequency log(s) 33 were run which may be converted in time 
domain, hiog(t). Using a nearby seismic trace heisft) 37 and the well log hiog(t) 33, a 
Combined Log and Seismic Response (CLSR) filter c(t) 39 can be obtained by solving 
Equation (4). Assuming that c(t) is invariant over the seismic survey, the high- 
resolution petrophysical property log response hj'(t) from every seismic trace, can be 
extrapolated (or predicted) 24 by solving Equation (5) to generate three-dimensional 
petrophysical property volume cubes (36 and 38) which have a much enhanced 
vertical resolution compared with the seismic data. 

[0047] The petrophysical property data produced by this invention may include 
but are not limited to the following data: Gamma Ray, V-shale, porosity, density, net- 
to-gross, resistivity, hydrocarbon saturation, permeability, and any combination. In 
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^eory. any peUophysical property ei*.r as a .wo^ata— log —ion, or a 
tt^dimensional ,og volume cnbe can be generated through the above proofs. 
,0048, Figure. 4(A) through 4(L) shows an experiment of the method that first 
erives a CLSR filter with a high-resolution log and then a low-r«olution s«sm,e 
^thensecondpr^iietsthegivenhigh-resolutionlogwiththeCLSR filter andlow- 

^«on seismie trace. Figures 4CA, and 4(B, are the high-.e<,uency og respons 
Snretion and it. « spec™., respectively. Figures 4(C, a.>d ^) a. * 
low-ftecuency se«mic wavelet and its ft^uency speC^s, respectrvely. F^es 
4(E, and 4(F, are hi^-fte,uency log fiom a weU and its fie,„ency . 
respectively. Figures 4(G) and 4(H) are the low-ftequency seismrc trace nea^ the wel 
I! ts fieUTspectrums, respectively. Figures 40) and 4(0 are the ^mbmed 
^^g-SeisnJc Response filter by solving B<,uafion (4) or (9), ^ - 
spitrums. Figures 4(K) and 4(L) are the predicted high—on log by ,oi™g 
^nation (5) or (10) and its fiouency spectrums. Note fi.t the ^cted^n^ 
llufion log and its spectrum in figures 4(K) and 4(L) are ahnost tdenfical to the 
original high-resolution log run ftom a well in figures 4 (E) and 4(F). 
,0049, AS shown in figure 4(H), the spect^m of the low-fiequency s^smic da« is 
ear zero when f.,uencies goes higher, and over which the high.fie,uency^g ^ 
in figure 4(0 has si^fican. amplimde in higher fie,uency. Thrs me^s that « . 
^e to do spectra, division as in the discussion about derivafion of E,uat.on (3)^ 
previous Jpossible tas. of spectral division is now solved using the — 

invlion for a stable log-seismic response filter c<., ™^ -thod provrdes a g^ 
predicfion of the high-resolution log as shown in fi^es 4(K) and 4(U, versus ^ 
original hi^-resolufion log run fiom a well in figures 4(E) and 4(F). V. pre^^ 
Z-luton log has enhanced resolution and can be used for 
glgical and geophysical analysi^mterpretaUon. geological modehng, and reservotr 

simulation. 



Express Mail No. ER888151472US 



